
function [Q,varargout] = alphaVarianceFull(yt,B,R,q,Q_0,T_0,prior)


Y = yt(:,q+1:end);
X =[];
for qx=q:-1:1    
    X = [X;yt(:,qx:end-q+qx-1)];
end

nu = R\(Y - B*X);
if size(yt,1)>53 || strcmp(prior,'diffuse')==1
[Q,sse] = drawVariance(nu,Q_0,T_0);
else
[Q,sse] = drawVarianceFull(nu,Q_0,T_0);
end

varargout(1) = {nu};
varargout(2) = {sse};


